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I.  INTRODUCTION 


With  a  continuing  interest  in  magnetic  anomaly  detec¬ 
tion  the  Naval  Postgraduate  School  has  studied  geomagnetic 
noise,  first  with  the  proton  precession  magnetometer  and 
later  with  the  optically  pumped  rubidium  type.  In  1972 
the  school  obtained  an  optically  pumped  cesium  magnetometer 
with  accessories  to  convert  it  to  a  gradiometer.  In  order 
to  take  full  advantage  of  the  inherent  high  frequency 
response  of  this  magnetometer  a  data  collection  system 
was  developed  to  observe  geomagnetic  and  gradient  noise. 
Figure  1.1  is  a  block  diagram  of  the  system. 

A  brief  introduction  to  the  earth's  magnetic  field 
and  other  measuring  instruments  is  presented  first  followed 
by  a  description  of  the  system  developed. 
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Figure  1.1  System  Block  Diagram 
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II.  THE  MEASUREMENT  FIELD 


A.  MAIN  FIELD 

World  magnetic  surveys  both  at  ground  level  and  more 
recently  by  satellite  indicate  that  the  principal  source 
of  the  main  geomagnetic  field  is  beneath  the  earth's  crust. 
Approximately  907®  of  the  main  field  exhibits  characteristics 
roughly  equivalent  to  the  field  of  a  terrestrially  centered 
short  bar  magnet  or  dipole  inclined  at  11.5°  to  the  axis  of 
rotation.  See  Figure  2.1. 


Figure  2.1  Dipole  Appearance  of  Geomagnetic  Field 
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The  field  strength  is  approximately  60  kilogammas  in 
the  polar  regions,  30  kilogammas  in  the  equatorial  region, 
and  50  kilogammas  in  the  neighborhood  of  37°  North  latitude 
on  the  west  coast  of  the  United  States.  Figure  2.2  shows 
a  world  chart  of  the  total  intensity  of  the  surface  field. ^ 
The  main  field  has  a  slow  secular  change  of  about  .1%  per 
year  and  it  is  a  regional  phenomenon  which  may  reach  .  37, 
in  certain  regions. 

B.  UNITS 

The  customary  unit  of  measurement  is  the  gamma,  regarded 
interchangeably  as  a  unit  of  induction  (gauss)  or  of  magnetic 
intensity  (oersted) . 

1  gamma  is  10~“*  gauss  (10-^  tesla) 
or 

1  gamma  is  10 oersted  (10  /4  ir  ampere  turn/meter) 

C .  ELEMENTS 

The  various  elements  of  the  magnetic  field  are  shown 
in  Figure  2.3.  They  are  defined  as: 

X  North  South  Component 

^Air  Force  Cambridge  Research  Laboratories,  Handbook  of 
Geophysics  and  Space  Environments,  p.  11-11,  1965. 
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Figure  2.2  World  Magnetic  Chart 


Y  East  West  Component 
Z  Vertical  Intensity 
D  Declination 

H  Intensity  of  the  Horizontal  Component 
F  Total  Intensity 
I  Inclination  or  Dip 

ZEIIITH  (-) 

..  t 

^  NORTH  (+) 


Figure  2.3  Elements  of  the  Field 


In  terms  of  the  more  familiar  bar  magnet:  A  magnet 
perfectly  free  to  turn  in  any  direction  would  show  the 
vector  direction  of  the  field.  It  is  almost  impossible 
to  suspend  a  magnet  like  that  so  two  are  used.  One  with 
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a  vertical  axis  (compass)  and  one  with  a  horizontal  axis 
(dip  needle) .  The  compass  shows  the  direction  of  the 
horizontal  component  H,  and  when  properly  oriented  the  dip 
needle  shows  the  inclination  or  dip. 

D.  TIME  VARIATIONS 

The  time  variations  of  the  geomagnetic  field  can  be 
broadly  categorized  into  quiet  variation  fields,  disturbed 
variation  fields ,  and  micropulsations . 

1 .  Quiet  Variation  Fields 

Geomagnetic  records  occasionally  show  a  very  smooth 
trace  that  clearly  indicates  a  pattern  of  daily  variation 
with  respect  to  solar  local  time.  This  pattern  is  named 
"solar  quiet  day  variation"  and  denoted  Sq.  There  is  also 
a  lunar  daily  variation,  denoted  L,  but  it  is  much  weaker 
and  more  variable.  Both  the  Sq  and  L  variations  result 
from  the  dynamo  action  of  the  atmosphere,  a  complex  phenom¬ 
enon  of  convective  motions  caused  by  the  solar  heating  of 
the  upper  atmosphere  under  the  influence  of  the  earth’s 
magnetic  field.  These  convective  motions  induce  ionospheric 
currents  and  electric  currents  in  the  earth  itself,  both 
of  which  cause  magnetic  variations.  The  effects  of  electric 
current  induction  within  the  earth  are  sometimes  referred 
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to  as  geological  disturbances  because  the  amount  of  induc¬ 
tion  is  dependent  on  the  composition  of  the  earth’s  crust. 
For  example,  the  intensity  of  induced  currents  in  the  ocean 
is  higher  than  those  in  the  continental  crustal  layers. 

Figure  2.4  shows  the  Sq  variation  recorded  in 
Monterey  during  a  quiet  period. 

2 .  Disturbed  Variation  Fields 

Any  fluctuation  other  than  the  quiet  day  solar  and 
lunar  variations  can  be  called  a  magnetic  disturbance. 

Other  than  disturbances  directly  attributable  to  a  solar 
flare  or  fluctuations  produced  by  upper  atmospheric  irreg¬ 
ular  motions,  all  are  caused  originally  by  disturbed  solar 
plasma  and  may  therefore  be  generally  classified  as  one 
family.  The  term  magnetic  storm  is  reserved  to  apply  to 
a  relatively  severe,  long  lasting  disturbance  with  recog¬ 
nizable  features,  but  all  other  irregular  variations 
deviating  from  the  Sq  and  L  variations  are  termed  simply 
magnetic  disturbances.  They  have  durations  lasting  from 
100  hours  to  several  minutes.  Figure  2.5  shows  a  magnetic 
disturbance  recorded  in  Monterey  that  immediately  followed 
the  Sq  variations  shovm  in  Figure  2.4. 
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Figure  2.5  Magnetic  Disturbances,  100  Gammas  Full  Scale 


3 .  Geomagnetic  Micropulsations 


Geomagnetic  micropulsations  are  field  variations 
with  periods  ranging  from  .2  seconds  to  approximately  10 
minutes.  A  terminology  recommended  by  the  International 
Association  of  Geomagnetism  and  Aeronomy  is  pc  for  pulsa¬ 
tion  continuous  and  pi  for  pulsation  irregular.  The  period 
range  for  these  are: 


pc 

1 

.2-5  sec 

pc 

2 

5-10  sec 

pc 

3 

10-45  sec 

pc 

4 

45-150  sec 

pc 

5 

150-600  sec 

Pi 

1 

1-40  sec 

Pi 

2 

40-150  sec 

A  general  picture  of  the  geomagnetic  field  in  the 
lower  frequencies  is  shown  in  Figure  2.6. 

Pc  1  pulsations  have  very  small  amplitude,  between 
.01  and  .1  gamma  and  have  a  sinusoidal  and  beating  form. 
The  pulsations  appear  in  bursts  and  usually  last  approxi¬ 
mately  30  minutes.  They  have  been  identified  by  their 
appearance  and  because  there  is  a  direct  relationship 
between  their  reoccurrence  period  and  mid-frequency.  No 
seasonal  or  latitude  dependence  has  been  identified  but 


21 


100  PERIOD1  (SEC)  01 


Figure  2.6  Geomagnetic  Micropulsation  Spectrum 

there  is  a  diurnal  variation  in  the  mid-frequency  events, 
high  frequency  events  occur  in  post  midnight  hours,  and 
low  frequency  events  in  the  early  afternoon.  Figure  2.7 
is  a  record  of  a  pc  1  micropulsation  observed  in  Monterey. 

Figure  2.8  shows  the  typical  form  of  the  pc  2,3 
pulsations.  They  also  exhibit  a  sinusoidal  and  beating 
character  similar  to  the  pc  1  type  but  their  amplitude  is 
slightly  greater,  usually  .5  gamma.  The  duration  of  the 
signals  seems  to  be  related  only  to  the  sensitivity  of  the 
detector,  which  has  lead  to  the  hypothesis  that  when  no 
signal  is  observed  the  pulsation  is  not  completely  absent 
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Figure  2.7  Pci  Micropulsation; 

1.8  gammas  full  scale;  45  second  record 
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but  rather  too  weak  to  be  detected.  Pc  2,3  events  occur 
mainly  during  daylight  hours  and  have  a  maximum  near  noon. 
The  average  period  undergoes  a  diurnal  variation,  being 
a  maximum  at  midday.  They  also  have  a  27  day  (solar  rota¬ 
tion)  periodicity  and  a  slight  seasonal  variation. 

As  seen  in  the  typical  pc  4  record  shown  in  Figure 
2.9,  these  micropulsations  slowly  rise  to  a  maximum  of 
between  5  and  20  gammas  and  then  symmetrically  decay.  The 
average  signal  lasts  about  one  hour.  In  mid-latitudes  a 
daytime  occurrence  pattern  can  be  observed  but  the  period 
is  shorter  in  the  daytime  than  at  night,  the  opposite  of 
the  pc  2,3. 

The  longest  period  variations,  pc  5,  display  ampli¬ 
tudes  from  10  to  hundreds  of  gammas .  Their  form  is  again 
sinusoidal  and  they  may  last  from  several  minutes  to  several 
hours.  The  largest  amplitudes  are  observed  in  the  auroral 
zone  and  the  occurrences  have  a  direct  relationship  with 
auroral  zone  phenomena. 

The  class  of  micropulsations  in  the  2  to  .005  Hz 
range  with  irregular  form  and  a  close  relationship  to  upper 
atmospheric  and  auroral  zone  phenomena  have  been  designated 
pi  1,2.  Their  amplitudes  vary  greatly,  the  larger  sometimes 
being  called  micropulsation  storms.  Figure  2.10  shows  the 
irregular  appearance. 
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Figure  2.9  Pc4  Micropulsation; 

7  gammas  full  scale;  6  minute  record 


Figure  2.10  Pi2  Micropulsation; 

7  gammas  full  scale;  6  minute  record 
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III.  MEASURING  INSTRUMENTS 


Today  a  wide  range  of  devices  are  available  for  measur¬ 
ing  the  magnetic  field,  ranging  from  the  extremely  old 
magnetized  needle  to  the  newer  superconducting  quantum 
types.  The  major  considerations  in  the  selection  of  a 
device  in  a  non-observatory  application  are  sensitivity, 
accuracy,  frequency  response,  mobility,  and  reliability. 
Three  types  of  magnetometers  that  have  developed  with  these 
considerations  as  guidelines  are  the  fluxgate,  optically 
pumped  resonance,  and  superconducting  quantum.  They  will 
be  described  and  a  comparison  of  the  different  varieties 
and  other  types  is  given  in  Table  I. 

A.  FLUXGATE 

A  fluxgate  magnetometer,  sometimes  referred  to  as  a 
saturable  core  or  saturable  inductor,  measures  magnetic 
fields  by  utilizing  the  variation  of  permeability  with 
intensity  of  magnetization  in  a  ferromagnetic  core.  An 
easily  saturable  ferromagnetic  core  of  high  permeability 
is  driven  cyclically  into  saturation  by  an  applied  magnetic 
field.  In  the  absence  of  an  external  field,  which  is 
usually  dc  or  very  low  ac  and  small  with  respect  to  the 
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applied  field,  the  induced  voltage  in  the  sense  windings 
will  be  symmetrical,  i.e.,  contain  only  odd  harmonics  of 
the  fundamental  driving  current  that  produced  the  applied 
field.  In  the  presence  of  an  external  field,  the  induced 
voltage  becomes  asymmetrical ,  the  asymmetry  and  the  corre¬ 
sponding  even  harmonics  being  proportional  to  the  amplitude 
of  the  component  of  the  external  field  parallel  to  the  axis 
of  the  core. 

The  fluxgate  magnetometer  has  been  used  widely  in 
various  applications  where  sensitivity,  stability,  and 
absolute  accuracy  has  been  sacrificed  for  the  features  of 
simplicity,  low  cost,  low  power,  and  ruggedness.  Since  it 
is  a  directional  device,  measuring  only  the  component 
parallel  to  the  axis  of  the  core,  a  triaxial  arrangement 
of  three  cores  is  necessary  for  total  field  measurements. 

B.  OPTICALLY  PUMPED  RESONANCE 

The  optically  pumped  resonance  magnetometers  (chiefly 
Rubidium,  Cesium,  and  Helium)  measure  the  total  intensity 
of  the  magnetic  field  by  an  extremely  accurate  detection  of 
the  Zeeman  effect.  Due  to  the  Zeeman  effect  the  energy 
levels  associated  with  a  valence  electron  become  split  into 
sublevels  whose  separations  are  dependent  on  the  total 


27 


intensity  of  the  magnetic  field.  The  detection  of  this 
splitting  is  based  on  the  principles  of  optical  pumping 
and  monitoring.  Optical  pumping  is  described  in  a  tutorial 
manner  in  Reference  4  and  for  the  Cesium  atom  in  Appendix  A. 

Figure  3.1  shows  a  simplified  block  diagram  of  a  typical 
pumping  apparatus.  A  light  source  producing  an  intense  col¬ 
limated  beam  of  light  is  first  filtered  and  circularly 
polarized  and  then  passed  through  a  resonance  cell  containing 
the  vapor  to  be  pumped.  As  the  light  is  absorbed  and  emitted 
by  the  sample  atoms,  optical  pumping  proceeds  and  results 
in  a  steady  state  non-equilibrium  population  distribution 
over  the  magnetic  states  of  one  of  the  ground  levels.  In 
this  pumped  condition  the  sample  is  prepared  for  resonance 
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Figure  3.1  Typical  Pumping  Apparatus 
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and  is  characterized  by  a  net  magnetization  M  =  Mz  along 
the  ambient  magnetic  field  direction.  Resonance  transitions 
between  the  Zeeman  states  of  the  pumped  level  are  induced 
by  subjecting  the  sample  to  an  RF  field  of  frequency  f-^ 
sufficiently  close  to  the  resonant  frequency  fQ  and  with 
component  H-^  perpendicular  to  the  ambient  field  Hq  =  Hz. 

The  RF  drive  field  introduces  Zeeman  coherence,  and  accord¬ 
ing  to  Bloch's  equations  the  steady  state  behavior  of  the 
magnetization  M  changes  and  is  no  longer  parallel  to  Hq  but 
now  has  a  transverse  component  Mj_  precessing  around  Hq  at 
the  frequency  of  H,  .  Consequently  the  Mz  component  along 
Hq  decreases,  and  the  light  transmission  also.  This  varia¬ 
tion  of  light  intensity  with  the  applied  RF,  being  a  maximum 
at  resonance,  can  be  used  as  a  detection  method,  or  the 
photocell  can  be  made  to  observe  the  signal  corresponding 
to  the  M  component,  which  is  precessing  or  rotating  at  the 
same  RF  applied  which  induced  the  resonance.  The  photocell 
output  is  then  at  the  resonant  frequency.  Appendix  B 
details  the  two  detection  systems  commonly  used.  The 
resonant  frequency  is  proportional  through  a  constant  to 
the  magnetic  field  (4.66  Hz/gamma  for  Rb;  3.499  Hz/gamma 
for  Cs;  28  Hz/gamma  for  He)  and  thus  the  field  measurement. 
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Figures  2.4  through  2.10  were  all  recorded  using  an 
optically  pumped  magnetometer  and  the  appearance  of  the 
micropulsations  differ  somewhat  from  those  one  would  find 
in  a  reference  handbook.  This  illustrates  an  undesirable 
feature  of  this  magnetometer  when  used  in  micropulsation 
work.  It  is  equally  sensitive  to  all  frequencies  and  the 
larger  amplitude,  lower  frequency  changes  such  as  the  di¬ 
urnal  variation  continually  drive  the  highly  sensitive 
instrument  off scale.  The  result  is  a  sloping  and  discon¬ 
tinuous  output  that  is  difficult  to  interpret. 

C.  CRYOGENIC  SQUID 

A  cryogenic  SQUID  (Superconducting  Quantum  Interference 
Device)  magnetometer  consists  of  three  basic  parts:  1)  a 
superconducting  ring  embodying  a  Josephson  junction;  2)  an 
externally  generated  HF  field  used  to  provide  a  small  field 
variation  superimposed  on  the  ambient  dc  field  and  therefore 
to  control  the  amount  of  flux  through  the  ring;  3)  an  elec¬ 
tronic  detection  system  to  observe  the  impedance  of  the 
coil  and  make  the  output  response  a  linear  function  of  the 
ambient  field. 

The  operation  utilizes  the  principle  that  if  the  flux 
through  the  ring  is  varied  on  either  side  of  a  minimum  or 
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maximum  energy  value,  associated  with  the  Josephson  coupling 
energy,  the  SQUID  will  generate  a  field  which  will  respec¬ 
tively  oppose  or  add  to  the  variation.  The  energy  contribu¬ 
tion  of  the  SQUID  aiding  or  opposing  the  applied  field 
variation  according  to  whether  the  ambient  dc  field  amounts 
to  a  minimum  or  maximum  energy  value  constitutes  a  para¬ 
metric  control  of  the  inductance  measured  at  the  coil 
terminals.  In  other  words,  the  SQUID  alternately  behaves 
diamagnetic ally  and  paramagnetically  and  will  decrease  and 
increase  the  inductive  reactance. 

Through  the  use  of  a  detector  to  monitor  the  coil  imped¬ 
ance  one  can  observe  the  variation  of  impedance  with  ambient 
magnetic  field.  Because  the  impedance  variation  is  periodic 
a  phase  locked  feedback  loop  is  normally  included  to  obtain 
a  response  linearly  related  to  the  field. 

Because  SQUID  compares  the  energy  values  of  the  ambient 
and  weak  superimposed  fields,  it  measures  only  a  change  in 
field,  not  the  total  intensity. 
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D .  GRAD  IOMETER 


A  gradiometer,  in  theory,  would  measure  the  space 
gradients  of  the  magnetic  field,  that  is: 

,  3x 

3x  3y  ,  3z 

3y_  ,  3y_  ,  3y 

3x  3y  3z 

dZ_  ,  dZ_  ,  3z 

3x  3y  3z 


Only  five  of  these  are  independent  because  X,  Y,  and  Z 
are  the  derivatives  of  a  potential  and  because  of  Laplace's 
equation.  In  practice  gradiometers  are  constructed  by 
taking  the  finite  difference  of  two  magnetometers,  there¬ 
fore  cancelling  the  ambient  field  that  is  approximately 
equal  at  both  sensors.  A  local  change  or  anomaly  that  is 
closer  to  one  of  the  two  sensors  is  not  cancelled  and 
therefore  recorded. 

The  disadvantage  of  the  gradiometer  is  that  with  respect 
to  a  single  sensor  magnetometer  most  magnetic  anomalies  fall 
off  as  the  inverse  cube,  1/r  where  r  is  the  distance  from 
the  anomaly  to  the  sensor,  but  the  gradiometer  anomaly  falls 
off  as  the  inverse  fourth  power,  1/r^  where  r  is  the 
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distance  from  the  anomaly  to  the  center  of  the  sensor 
configuration. 

In  Reference  10,  Langan,  Lynch,  and  Slack  reported  a 
gradiometer  used  in  geophysical  prospecting  that  measured 
the  vertical  gradient  A  H/az  where  H  was  the  magnitude  of 
the  anomaly  and  z  the  distance  to  the  anomaly.  It  was 
constructed  of  two  optically  pumped  sensors  vertically 
separated  and  towed  by  a  helicopter. 

Both  the  optically  pumped  and  SQUID  magnetometers  are 
easily  converted  to  gradiometers  because  of  the  form  of 
their  output,  but  they  will  measure  different  gradients. 
The  optically  pumped  gradiometer,  constructed  of  total 
field  magnetometer  sensors  will  measure  the  total  field 
gradient  AF/  Az.  A  SQUID  gradiometer,  constructed  of 
devices  sensitive  only  to  a  field  change  will  measure  a 
field  component  gradient  AZ/  Az,  for  example. 
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Table  1  Comparison  of  Magnetometers 


IV.  SYSTEM  DESCRIPTION 


The  overall  system  can  be  broken  into  three  subsystems: 

1)  the  data  collection  subsystem;  2)  the  periodmeter  and 
interface  subsystem;  3)  the  processing  and  display  subsystem. 


A.  DATA  COLLECTION  SUBSYSTEM 

The  magnetometer  used  was  the  Varian  V-4938  Cesium 
system.  It  used  an  optically  pumped  cesium  sensor  of  the 
auto-oscillating  type.  A  block  diagram  of  the  Varian  sys¬ 
tem  is  shown  in  Figure  4.1. 


SENSOR 


Figure  4.1  Varian  Magnetometer 
The  sensor  head  produces  the  RF  signal  (Larmor  frequency) 
directly  proportional  to  the  total  magnetic  field  intensity. 
This  signal  together  with  the  input  power  for  the  sensor  is 
conducted  along  a  200  foot  coaxial  cable  to  the  main  readout 
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unit.  The  readout  unit  converts  the  Larmor  frequency  to  a 
highly  stable  dc  voltage  which  drives  the  pen  recorder. 

The  system  was  located  in  a  wooded  area  approximately  one 
mile  from  the  Naval  Postgraduate  School,  which  provided  a 
near  noise  free  environment. 

A  more  complete  block  diagram  of  the  V-4938  system  is 
shown  in  Figure  4.2.  The  RF  Larmor  frequency,  which  was 
approximately  176  kHz,  from  the  sensor  is  mixed  in  the  read¬ 
out  unit  to  the  100  Hz  to  3  kHz  range  and  then  applied  to 
a  highly  stable  frequency  to  voltage  converter.  This  dc 
voltage,  which  now  represents  the  field  value,  is  passed 
through  biasing  circuitry  to  interface  it  to  the  pen  recorder. 

The  system  configuration  used  for  data  collection  differed 
from  the  Varian  recommended  system  of  Figure  4.2.  This  was 
necessary  because  of  the  requirement  to  record  the  output 
on  magnetic  tape.  The  actual  system  used  is  shown  in  Figure 
4.3.  This  configuration  was  used  because  the  mixer  output 
of  the  readout  unit  was  a  square  wave  and  its  frequency 
content  exceeded  the  5  KHz  range  of  the  tape  recorder.  The 
coupler  unit  shown  connects  between  the  readout  unit  and 
the  sensor  to  provide  dc  power  to  the  sensor  and  it  also 
contains  a  mixer,  but  its  output  is  a  sinusoid  rather  than 
a  square  wave  and  recorded  nicely.  The  coupler  units  primary 
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Figure  4.2  Magnetometer  Signal  Flow 
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Figure  4.3  Magnetometer  Data  Collection  Configuration 


purpose  is  to  convert  the  system  to  a  gradiometer  which  is 
shown  in  Figure  4.4. 


Figure  4.4  Gradiometer  Data  Collection  Configuration 

Again  this  is  not  the  configuration  recommended  but 
this  provided  both  a  magnetometer  and  gradiometer  output 
suitable  for  recording  and  still  allowed  the  pen  recorder 
to  monitor  the  gradiometer  signal.  The  bias  applied  to 
the  one  sensor  was  required  to  generate  a  difference 
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frequency  between  the  sensors  sufficient  to  drive  the  mixer 
circuitry.  The  sensors  were  located  approximately  180  ft 
apart  resulting  in  a  Larmor  frequency  difference  of  only 
approximately  30  Hz,  and  the  mixers  required  a  difference 
of  at  least  100  Hz  and  optimally  between  1  and  3  kHz. 

Figure  4.5  shows  the  readout  units  mixer  output  with  an 
insufficient  difference  frequency  and  Figure  4.6  shows  the 
coupler  units  sinusoidal  mixer  output. 

One  other  system  characteristic  imposed  a  limitation, 
that  was  the  orientation  sensitivity  of  the  single  cell 
sensor.  For  optimum  performance  the  sensor  head  needed  to 
be  oriented  both  toward  magnetic  North  and  45°  to  the  total 
field  vector.  Figure  4.7  shows  the  non -metallic  housing 
and  Figure  4.8  the  geometry. 

Figure  4.9  shows  the  system  components  and  equipment 
used  to  monitor  and  record  the  data. 
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Figure  4.5  Distorted  Readout  Unit  Mixer  Output 


Figure  4.6  Coupler  Unit  Mixer  Output 
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Figure  4.7 


Non -Metallic  Housing  Used  to  Fulfill 
Orientation  Requirements 
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field 


Figure  4.8  Sensor  Housing  Geometry 
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B.  PERIODMETER  AND  INTERFACE  SUBSYSTEM 


The  frequency  response  of  the  magnetometer  when  used 
with  the  readout  unit  is  limited  to  20  Hz  by  the  discrimin¬ 
ator  and  even  more  by  the  pen  recorder  response.  The 
theoretical  frequency  response  of  the  sensor  is  limited 
only  by  the  bandwidths  of  the  amplifier  and  feedback  loop, 
which  are  extremely  wide.  To  take  advantage  of  this  the 
mixer  output  was  recorded  on  an  HP-3960  FM  tape  recorder. 

This  signal  still  contained  the  high  frequency  components 
present  in  the  Larmor  signal.  A  periodmeter  and  interface 
unit  were  constructed  to  measure  the  period  of  successive 
cycles  and  interface  this  digital  count  into  the  XBS-9300 
computer,  for  further  display  on  the  AGT-10  graphics  console. 

Figure  4.10  is  a  block  diagram  of  the  periodmeter.  The 
data  signal  is  passed  through  a  simple  RC  high  pass  filter 
to  reject  any  60  Hz  modulation  that  might  be  present  and 
then  input  to  a  zero  crossing  detector.  The  falling  edge 
of  this  is  used  to  trigger  a  one-shot  that  inhibits  the 
otherwise  free  running  4  MHz  count.  The  one-shot  pulse 
is  inverted  and  delayed  and  the  rising  edge  used  to  trigger 
flip  flop  type  memories.  Once  the  count  is  latched  another 
one  shot  is  used  to  clear  the  count  and  provide  a  signal 
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Fig  4.10  PERIODMETER 


indicating  the  cycles  completion.  The  cycle  then  repeats. 

The  actual  circuits  and  timing  diagrams  are  presented  in 
Appendix  C. 

A  direct  access  channel  into  the  XDS-9300  was  used  to 
input  the  binary  count  directly  into  sequential  memory  loca¬ 
tions.  Figure  4.11  is  a  block  diagram  of  the  periodmeter 
and  interface.  To  implement  the  input  the  following  sequence 
of  instructions  were  used: 

EOM  032400 
SKS  032400 
BRU  $-1 
PIN  ITEM 

With  the  execution  of  the  EOM  instruction  the  contents  of 
the  control  register  are  available  at  the  output  connector. 
This  is  decoded  and  toggles  a  flip  flop  to  alert  the  period- 
meter  that  the  XDS  is  ready  to  accept  an  input.  This  EOM 
alert  and  the  periodmeter ' s  completion  signal  are  logically 
combined  and  send  a  ready  signal  to  the  XDS-9300.  Upon 
receipt  of  this  ready  signal  and  the  execution  of  the  PIN 
the  data  is  transferred  and  stored  in  the  location  specified 
in  the  PIN  instruction.  The  SKS  instruction  is  used  to  sense 
if  both  the  computer  and  periodmeter  are  ready.  If  both  the 
EOM  alert  and  completion  signal  are  not  present,  the  next 
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Figure  4.11  Periodmeter  and  Interface  Block  Diagram 


instruction,  branch  unconditionally  back  to  the  SKS ,  is 
executed.  This  looping  continues  until  both  systems  are 
ready.  If  both  are  ready  the  branch  instruction  is  skipped. 
The  ready  signal  is  also  used  to  gate  the  data  through  in¬ 
verters  because  it  is  inverted  prior  to  being  transferred 
into  the  XDS-9300.  The  actual  circuit  diagrams  are  pre¬ 
sented  in  Appendix  C. 


C.  DISPLAY  SUBSYSTEM 

A  block  diagram  of  the  data  flow  in  the  interface  and 
display  subsystems  is  shown  in  Figure  4.12. 
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Figure  4.12  Data  Flow  From  Periodmeter  Through  Display 


The  physical  block  length  of  data  input  to  the  XDS-9300 
was  12,288  points.  A  mean  and  standard  deviation  of  these 
points  was  found  and  values  not  within  five  standard 
deviations  were  discarded.  This  eliminated  erroneous  data 
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points  caused  by  noise  spiket  either  in  the  computer, 
periodmeter,  or  on  the  recording.  The  block  was  then 
buffered  out  in  groups  of  4096  points  and  written  on 
magnetic  tape.  To  be  displayed  the  data  was  buffered 
back  in,  again  in  groups  of  4096  points,  normalized,  and 
displayed  on  the  AGT-10  graphics  console.  Both  programs 
are  self-explanatory  and  are  presented  in  Appendix  D. 
Signal  flow  diagrams  for  the  two  are  shown  in  Figure  4.13 
and  4.14. 
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Figure  4.13  Signal  Flow  Diagram  of  Input  Program 
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Figure  4.14  Signal  Flow  Diagram  of  Display  Program 
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V.  RESULTS  AND  CONCLUSIONS 


The  system  was  tested  in  three  phases:  1)  laboratory 
signals;  2)  tape  recorded  laboratory  signals;  3)  magnetom¬ 
eter  signals. 

First,  constant  frequency  signals  from  a  laboratory 
generator  were  applied.  The  frequencies  sampled  and  the 
resulting  period  measurement  are  listed  in  Table  2. 

FREQUENCY  Hz  MEAN  PERIOD  ^sec  THEORETICAL  PERIOD  ^sec 


1000.0 

999.25 

1000.00 

1000.1 

999.00 

999.00 

1000.5 

998.75 

999.50 

1001.0 

998.25 

999.00 

1003.0 

996.25 

997.01 

1005.0 

994.25 

995.02 

Table  2 


The  .75  jj  sec  difference  between  theoretical  and  actual  is 
due  to  the  clear  and  latch  time  of  the  periodmeter.  These 
measurements  were  within  +1  count  or  +25  ^sec  with  a  4MHz 
time  base.  This  established  the  resolution  of  the  device 
as  . 1  Hz  or  better.  The  outputs  displayed  for  the  sampled 
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frequencies  of  1000.0  Hz,  1000.1  Hz,  and  1001.0  Hz  are 
shown  in  Figures  5.1  -  5.3  and  a  distinguishable  difference 
is  noticeable.  To  insure  the  system  would  follow  a  rapidly 
changing  frequency  a  swept  frequency  signal  was  applied  and 
the  display  shown  in  Figure  5.4. 

The  second  phase  of  testing  was  designed  to  determine 
the  magnitude  of  the  tape  recorder  flutter  and  noise.  It 
involved  recording  the  same  frequencies  used  in  phase  1  and 
comparing  the  results.  The  frequencies,  mean  period  measure¬ 
ment,  and  standard  deviation  for  each  are  listed  in  Table  3. 

FREQUENCY  Hz  MEAN  PERIOD  usee  STANDARD  DEVIATION  ^sec 


1000.0 

999.25 

.575 

1000.1 

999.00 

.700 

1000.5 

998.75 

.850 

1001.0 

998.25 

.550 

1003.0 

996.25 

.650 

1005.0 

994.25 

.675 

The  mean  period,  calculated  for  12,288  periods,  remained 
the  same.  The  average  standard  deviation  was  .667  jjsec. 
This  established  the  RMS  noise  level  of  the  system  when 
operated  in  the  1200  Hz  range  to  be  approximately  1  Hz  or 
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Figure  5.1 


Output  Display  of  Direct  Input 
1000.0  Hz  Laboratory  Signal 


1000. 


PERIOD 
p  SEC 


Figure  5.2  Output  Display  of  Direct  Input 
1000.1  Hz  Laboratory  Signal 
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1000 


F igure  5 . 3 


Output  Display  of  Direct  Input 
1001.0  Hz  Laboratory  Signal 


TIME,  SEC 


1000. 


PERIOD 

juSEC 


990. 

0. 


Figure  5.4  Output  Display  of  Direct  Swept  Frequency 
Laboratory  Signal 


56 


.27  gamma.  The  AGT-10  displays  for  the  sampled  frequencies 
of  1000.0  Hz  and  1000.1  Hz  are  shown  in  Figures  5.5  and  5.6. 
The  previously  discernible  difference  is  now  buried  in  tape 
recorder  noise.  Figure  5.7  shows  the  frequency  transition 
from  1001.0  to  1003.0  Hz  and  this  2  Hz  difference  is  easily 
distinguishable. 

The  final  phase  of  testing  involved  introducing  a  known 
magnetic  disturbance  and  analyzing  the  magnetometer  output. 

A  72  turn  Helmholtz  coil,  in  the  configuration  of  Figure  5.8, 
was  used  to  introduce  the  disturbance.  The  sensor  was 
placed  on  the  axis  of  the  coil  and  the  coil  aligned  with 
the  total  field  vector.  This  insured  the  induced  field  of 
the  coil  would  directly  aid  or  oppose  the  ambient  field 
and  allowed  changes  of  predetermined  amounts.  With  the 
simplified  geometry 


H  = 


N  I  R" 


4  7T 


10 


-Z 


where  N  =  number  of  turns  in  the  coil 


I  =  current  in  the  coil 


R  =  radius  of  the  coil 


r  =  distance  from  coil  to  sensor 


was  required  for  conversion  from 
ampere-turns  per  meter  to  gammas. 
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990 
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Figure  5.5 


Output  Display  of  Recorded 
1000.0  Hz  Laboratory  Signal 
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990. 


Figure  5.6 


Output  Display  of  Recorded 
1000.1  Hz  Laboratory  Signal 
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TIME, SEC 

Figure  5.7  Output  Display  of  RecordedlOOl. 

and  1003.0  Hz  Laboratory  Signal 
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MAGNETIC 

NORTH 


Figure  5.8  Coil  and  Sensor  Geometry  Used 
to  Disturb  the  Field 
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The  coil  utilized  had  a  mean  radius  of  33  cm  and  was  placed 
3.5  m  from  the  sensor.  Table  4  tabulates  the  various 


currents  and  observed  versus  predicted  changes.  Figures 
5.9  through  5.30  show  the  magnetometer  pen  recorder  output 


and  a 

selected  time 

frame  of  the 

corresponding 

AGT-10  display 

THEORETICAL 

MAGNET. 

THEORETICAL 

DISPLAY 

I  ma 

AH 

AH 

Af 

Af 

400 

45.964 

46.0 

160.874 

160.893 

300 

34.473 

35.0 

120.655 

120.666 

240 

27.578 

27.5 

96.524 

96.531 

200 

22.982 

22.7 

80.437 

80.442 

160 

18.386 

18.7 

64.350 

64.354 

120 

13.789 

13.8 

48.262 

48.264 

100 

11.491 

11.1 

40.218 

40.219 

90 

10.342 

9.6 

36.197 

36.198 

80 

9.193 

9.0 

32.175 

31.176 

70 

8.044 

8.4 

28.153 

28.154 

60 

6.895 

6.8 

24.131 

24.132 

40 

4.596 

4.2 

16.087 

16.086 

20 

2.298 

2.1 

8.044 

8.045 

10 

1.149 

1.1 

4.022 

4.023 

5 

.574 

.65 

2.011 

2.012 

Table  4 
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FIELD  CHANGE 


Figure  5.9  Pen  Recorder  Output  for  AH=  45  gamma 


TIME,  SEC 

Figure  5.10  AGT-10  Display  for  AH=  45  gamma 
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FIELD  CHANGE 


Figure  5.11  Pen  Recorder  Output  for  AH=  35  gamma 
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AGT-10  Display  for  A H=35  gamma 
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Figure  5.12 
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FIELD  CHANGE 
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Figure  5.13  Pen  Recorder  Output  for  A  H=  28  gamma 
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Figure  5.14 


AGT-10  Display  for 


AH=  28  gamma 
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FIELD  CHANGE 


Figure  5.15  Pen  Recorder  Output  for  aH=  18  gamma 
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Figure  5.16 


AGT-10  Display  for 


A  H=  18  gamma 
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FIELD  CHANGE.  C'u- 


Figure  5.17  Pen  Recorder  Output  for  AH=  11  gamma 
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Figure  5.18  AGT-10  Display  for  aH=  11  gamma 
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FIELD  CHANGE 


TIME,  SEC 

Figure  5.19  Pen  Recorder  Output  for  A  H=  8  gamma 
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0.  1.2 
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Figure  5.20  AGT-10  Display  for  A H=  8  gamma 
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FIELD  CHANGE 


Figure  5.21  Pen  Recorder  Output  for  A  H=  4  gamma 
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Figure  5.22  AGT-10  Display  for  AH=  4  gamma 
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FIELD  CHANGE 


Figure  5.23  Pen  Recorder  Output  for  AH=  2  gamma 
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Figure  5.24  AGT-10  Display  for  AH=  2  gamma 
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FIELD  CHANGE 


Figure  5.25  Pen  Recorder  Output  for  AH=  1  gamma 
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Figure  5.26  AGT  Display  for  AH=  1  gamma 
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FIELD  CHANGE 
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Figure  5.27  Pen  Recorder  Output  for  AH=  .6  gamma 


1000. 


0.  1.2 

TIME,  SEC 

Figure  5.28  AGT-10  Display  for  AH=  .6  gamma 
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To  obtain  a  precise  measurement  of  the  magnetometer’s 
response  time  the  AGT-10  display  was  expanded.  Figure  5.29 
shows  the  rising  edge  and  Figure  5.30  the  falling  edge.  As 
seen  the  rising  edge  required  approximately  three  periods 
or  2.5  msec  to  respond  and  the  falling  edge  only  one  period 
or  .8  msec. 


TIME,  SEC 

Figure  5.29  Expanded  View  of  Rising  Edge 
of  Magnetometer  Response 
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TIME,  SEC 

Figure  5.30  Expanded  View  of  Falling  Edge 
of  Magnetometer  Response 
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The  periodmeter  and  interface  unit  measured  the  period 
of  the  signal  with  an  absolute  accuracy  of  .75  jj.  sec  and 
properly  input  the  value  to  the  computer.  The  resolution 
of  the  device  was  .25  jj,  sec  using  a  4  MHz  time  base  and 
this  could  be  increased  by  using  a  faster  time  base.  The 
system  was  successfully  operated  with  a  10  MHz  clock  fre¬ 
quency  and  resolution  of  .1  j^sec.  The  use  of  an  increased 
clock  frequency  or  averaging  of  periods  to  increase  the 
periodmeters  resolution  would  not  improve  the  system  sensi¬ 
tivity.  The  overall  system  sensitivity  is  limited  by  the 
flutter  and  noise  of  the  tape  recorder.  Two  methods  were 
attempted  to  overcome  this:  1)  the  use  of  a  synchronized 
oscillator  constructed  of  two  multivibrators  synchronized 
with  the  zero  crossings  of  the  input  signal;  2)  recording 
a  1  kHz  reference  signal  on  a  second  channel  of  the 
recorder  and  multiplying  it  up  for  use  as  a  time  base. 

The  pulse  width  instability  of  the  multivibrators  was 
greater  than  the  flutter  noise  so  the  synchronized  oscillator 
was  unsatisfactory.  A  phase  lock  loop  type  frequency  multi¬ 
plier  was  constructed  to  multiply  the  1  kHz  to  1  MHz  but 
the  delay  time  through  the  loop  was  too  great  to  cancel  the 
flutter  for  a  period  to  period  measurement. 

To  achieve  greater  system  sensitivity  and  resolution  some 
means  of  flutter  cancellation  is  necessary. 
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APPENDIX  A 


OPTICAL  PUMPING  OF  THE  CESIUM  ATOM 


Figure  A-l  shows  the  energy  levels  of  the  cesium  atom. 

Absorption  of  photons  of  light  by  the  cesium  valence  elec- 

o 

trons  cause  their  energy  to  be  increased  from  the  ground 

2 

state  to  the  Pi  state.  Due  to  the  light's  polarization 

the  magnitude  of  the  quantum  number  m  is  also  increased  by 

one.  An  excited  atom  quickly  emits  its  energy  and  returns 

to  any  of  the  ground  state  sublevels.  The  absorption  and 

re-emission  is  repeated  resulting  in  a  shift  to  the  m=+4 

2 

sublevel.  Since  the  Pi  state  contains  no  m  level  above 

h 

2 

+4,  transitions  out  of  the  Si ,  m=+4  level  are  not  possible. 

*2 

This  is  precisely  optical  pumping  and  under  these  conditions 
there  can  be  no  further  absorption  of  light. 
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Figure  A-l  Energy  Levels  of  Cesium  Atom 
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APPENDIX  B 


OPTICAL  DETECTION  SYSTEMS 


The  first  method  used  to  obtain  a  continuous  record  of 
the  total  field  intensity  used  a  control  loop  arrangement 
and  was  therefore  called  a  Controlled  Magnetometer.  A 
block  diagram  of  it  is  shown  in  Figure  B-l. 


DIRECTION  OF  FIELD 

^ - ► 


FILTER  RESONANCE 

CELL  PHOTOCELL 


A  small  frequency  modulation  is  applied  to  the  local  RF 
oscillator  to  sweep  it  back  and  forth  through  resonance. 
Owing  to  this  sweeping  the  signal  collected  at  the  photo¬ 
cell  is  also  modulated  at  the  same  frequency,  but  its  phase 
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with  respect  to  the  initial  modulating  signal  depends  on 
the  position  of  the  local  RF  with  respect  to  the  resonant 
frequency.  The  photocell  output  and  modulating  signal  are 
phase  detected  and  the  output  is  used  as  an  error  signal 
in  a  control  loop  to  maintain  the  local  RF  oscillator  coinci- 
cent  with  the  resonant  frequency.  This  system  does  not  place 
severe  restrictions  on  the  performance  of  the  photocell  or 
its  amplifier,  has  high  accuracy,  good  stability,  and  good 
low  field  capability. 

One  attribute  the  controlled  system  does  not  have  is  a 
good  response  to  transients.  Figure  B-2  shows  the  block 
diagram  of  a  system  that  has  the  excellent  response  to 
transients  necessary  to  detect  and  observe  micropulsations. 


DIRECTION  OF 
FIELD 


Figure  B-2  Self -Oscillating  Magnetometer 
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In  theory  a  second  light  beam  is  placed  at  right  angle  with 
respect  to  the  first  and  the  variation  of  the  M  component 

A 

of  the  magnetization  is  observed.  Since  this  transverse 
component  oscillates  at  the  frequency  of  the  signal  injected 
and  the  amplitude  is  a  maximum  at  resonance,  the  photocell 
is  simply  amplified,  phase  shifted,  and  re-injected  to  the 
coils.  The  system  will  then  oscillate  at  the  resonant 
frequency.  In  practice  it  is  possible  to  use  a  single  light 
beam  inclined  at  45°  to  the  direction  of  the  field  to  be 
measured.  This  places  an  orientation  restriction  on  the 
sensor,  but  systems  employing  two  and  three  resonance  cells 
have  been  developed  that  overcome  this.  This  system  places 
severe  restrictions  on  the  amplifier  and  especially  on  the 
accuracy  of  the  phase  shift,  which  is  required  to  maintain 
oscillation.  It  has  a  low  absolute  accuracy  and  poor  long 
term  stability  but  as  previously  mentioned,  excellent 
response  to  transients. 


79 


7407 


APPENDIX  C 


PERIODMETER  AND  INTERFACE  CIRCUITS 
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Figure  C  l  Counting  Circuitry 
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Figure  C-2  Input  Circuitry  Timing  Diagram 
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Figure  C3  Input  Circuitry 
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FIGURE  C4  TIME  BASE 
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FIGURE  C  5  COMPUTER  INTERFACE 
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APPENDIX  D  COMPUTER  PROGRAMS 


-joe 

->FCFTP  A  N  LS,GC,S 

DIMENSION  CGUNT  ( 1 3000  )  » I  VAR (6  ) 
INTEGER  CO-NT 
CA7A  N$AMP,MULT/4G96,3/ 
L£LK=NSAM PENULT 
106  FCRNAT (  IX, 12 ( 110)  ) 

999  CONTINUE 

INITIALIZE  VARIALBLSS 

I  N  /S  >=  1 
I  MN  =  100000 
SLMO. 

S IN  SC  =0  • 

C  =  C. 

I E  A  C  =0 

eAc=c. 


INPUT  CATA  FROM  PE RICCNETER/ INTERFACE 


CC  1  =  1, LI 

EON  C32400 
SKS  032400 
EFL  $-1 
FIN  ITEN 
ICCLNT  =0 

I C  C  LNT=  LA  N  C ( ITE^ , 177B) 

I CCLNT=  L ICR ( I  COUNT ,LANC ( LRS ( ITEN » 1  )  ,37600  6  )  ) 
ICCLNT=LICR(ICC'UNT,  LANC(LRS(  I  TEN  ,  1)  ,140000  6  )  ) 

CCLNT (I  )  =  I  COUNT 

2  CONTINUE 

FCFN  AN  INITIAL  AND  THEN  PR6CICE  MEAN  £  SDEV 

CC  3  I  =  1 »  L6LK 
SIMSUK  +  CCLNT  (I  ) 

SLNSC=SLMSC+ COUNT  (I  )**2 
C  =  C-H.O 

3  CONTINUE 

M  E  AN  =  SU  N/C 

SCEV=SQRT (A6S (C*SU*SG-SLM**2 )  )/C 
V  P ITE (6 ,100 )MEAN,SCEV 
ICC  FCFNAT( '0  MEAN  = ' , F 12. 4/ 

1*  SCEV=',F12.4) 

00  4  I  =  1 , LBLK 

1F(ABS(CCUNT( I)-MEAN} ,LE.5.*SCEV)  GC  TO  4 
SIMSUM-CCLNT  (I  ) 

S INSQ=S LNSC-COUNT (I )**2 
E  AC=BAC  +1  • 

4  CONTINUE 
C=C-6AC 

N  E  A N=  SUN/C 

SCEV=SQFT  <AES(C*SUf'SC-SLM**2)  )/C 
W  F ITE (6 *  10 1 ) HE AN  * SCEV 

101  F  CFNAT ( *0  TRUE  M£AN  = » , F12.4/ 

1'  TRUE  SDEV=* ,F12.4  ) 

FOR  N  A  NEW  ARRAY  CONSISTING  OF  POINTS  WITHIN  10  SDEV 

I E  A  0=0 
N  £  W I =0 

CC  6  1  =  1, LBLK 

IF (ABS(CCLNT(I ) -MEAN ) .L E. 10.*S DE V ) NEW  =  NE W I  + 1 ;  GC  TO  7 
I EAC=I6 AC+1 

WRITE(6,1C2)  18 AD , I , CCLNT ( I  ) 

102  FCFNA.7(«0  THIS  IS  BAC  PCINT  NO‘,15/ 

1’  ITS  OLD  SUBSCRIPT  WAS', 18/ 

2*  ITS  VALUE  WAS • , 1 10 ) 

C-C  TC  6 

7  CCLNT(NEWI  )=C0UNT( I  ) 
t  CONTINUE 


I = 1 , LBLK 
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OOO  OOO  OOO  ooo 


CC  15  I =NEWI+1 > LBLK 
CCUTII  )  =  G 
15  CONTINUE 

WRITE (6  ,103)  NEW! 

102  FC  F  N AT ( 'C  NEW  1=' ,110) 

ZERO  FILL  AFRAY  FROM  LAST  POINT  TO  ENC  OF  BLOCK 

CC  6  I  =  NEW  1  +  1 , LBLK 
CCLNTII  )  =0 
E  CCNTINLE 

FINC  A  MINIMUM  AND  MAXIMUM  VALLE 

CC  9  I = 1, N  EW  I 
ININ=AMINQ<CCUNT< I ) , I  MINI 
IMAX=AMAXO<CGLNT( I ) , IMAX5 
9  CCNTINLE 

W  R I TE  ( 6 » 104  )  IMI N  » ! M A  X 
104  FCFMATt'G  MINIMUM* *, 1 10/ 

1*  MAXIMUM*' ,110) 

ELFFEFCUT  TC  TAPE  THE  ARRAY  IVAR  CONTAINING  VARIABLES 

I V A  R ( 1 ) =NEW I 
I\AR(2)=MEAN 
I  V A  P ( 3 ) *  SCE V 
IVAR{4)*IMIN 
IVAR<5) =IMAX 

WRITE ( 6 ,105 ) ( IVAR ( I  )  ,1=1,5) 

1 C 5  FORMAT!  1X.5UO) 

CL7FLTI 102) 'WRITE  ON  TAPE* 

REAC( 1C1,601)IANSS 
601  FCFMAT ( A4 ) 

IF(  IANSS.NE.3HY5S)  C-0  TO  603 
CALL  eUFFERCUTt 1, 1 « I  VAR ,5 , 1 NC ) 

11  IF(INCoEC.l)  GO  TO  11 

ELFFERCUT  TC  TAPE  THE  ARRAY  CGLNT  IN  GROUPS  OF  4096 

CC  12  1*1, MULT 
J  =  1-1 

K  =  ( J=*NS  AM  F  )  +  1 

CALL  6UFFERCUT(1,1,CCUNT(K)  ,NSAMF,IND) 

12  1  F (  INC. EC. 1)  GO  TC  13 
12  CCNTINLE 

6 C 2  CCNTINLE 

CLTFUT  1 102  )  • REAO  ANOTHER  14336' 

REACdCl,  600  )1  ANS 
6CC  F CF  MAT ( A  4 ) 

IF  ( IANS  .EC.3FYES  )  GO  TO  999 
ENC 

^LCAC  X  F , MA  P 
->CA  T  A 


86 


onn  nnn 


SPA7CI- 
1>>  C  A7A 
C07C67 
$  EN  C 
-AG7 
->  JC  E 

-FORTRAN 


CC1C6711 


flbriflch) 

OUTPUT ( 10 1 ) I  SR , • FI* 
GG  TC  2 


AN  LS.GC 

DIMENSION  I  AX  IS (10), INCAT  A (4105) ,  DATA (4 10 1 
11  C  C  I R ( 5  ), ITCIR< 15 ) » ITEXT (24) ,  T (2 ) 

INTEGER  R 5 S P ( 4 ) 

CATA  NSAMP,MULT,NCIS/4G96,3, 16/ 

I  C  E  V  =  1 

LELK=NSAMP*MLL7 
CIV*FLGA7 (NCIS) 

N  EL  K  =  1 

I  AXIS  (-1  )  =  I F  E  A  D  (  0 ,  10  ) 

I  AX  IS  <  2 )= IPACK (-• 5  » .75 ,0 ) 

I  A  X  IS (3 )= I P ACK ( -.5 , 1 • ,0  ) 

I  A  X  IS  (  4  ) = I  PACK ( - • 5  »C • , 1 ) 

I  AXIS (5 >  = IPACK(0. ,C., 1 ) 

I  AXIS (6 >  =  IPACK <0.  ,-.01,1) 

I A  X  IS ( 7 )  =  I  PACK ( 0 • , 0. , 0  ) 

IAXIS18 )  =  IPACK( .5,0.,  I  ) 
IAXIS(9)=iPACK(.5,-.01,l) 

999  CENT  INC E 
T IF  E=C. 

LFEC=NSAMP 

2  CALL  FN  S ( ICEV » ISW » IER  ) 

If  (  IER.NE.O 
IFUSW.NE.Q) 

BLFFEF IN  ARRAV  IVAR  KITH  VARIABLES 

CALL  EUFFEPIN(1,1,IVAR, 5, INC) 

2  1  F  (  IN  C.  EG  .  1 )  GO  TC  3 
NEfcl=lVAP.(l) 

NEAN  =  FLCAT( IVAR(2 )  )/4. 

S  CEV=FLCAT ( IV Aft (2  )  )/4. 

N IN  =  FLQAT (  IVAR(4)  )/4. 

MAX=FLOAT ( IVAR( 5 ) )/4. 

WFITE(6»8CG)( IVAR ( I  )  ,1  =  1,5) 

800  FCFNATllX, 'NO.  CF  POINTS* ', I 10/ 

1*  N  EAN  = ',110/ 

2  *  SEEV=',I10/ 

2'  MNIMLM**  ,  110/ 

4'  NAXIMLM=  '  , 110) 

V  R I T  E  M  S  AN , S  CcV  , MIN, VAX  CN  SCREEN 

CALL  CGIN  IT  t IDEV,  IGCIR, 5, IER) 

IF(  lER.NE.OJOLTPUT  (  10  1  )IER,  'CGI' 

CALL  0TINIT(IDEV,ITCIR,15,IER) 

IF ( IER.NE.O  OUTPLT ( 101 )  IER,'D2' 
ENCCDE(96, 200, ITcXT) 

FCF  NAT (  ' PE° IGD' ) 

CALL  TEXTC(IDEV, ITcXT ,24,3,26, 1,2, IER) 

IF( IER.NE.O ) OUTPLT < 101) IER, 'PER  ICC' 
ENCCCE(96, 201, ITEXT) 

201  FC  F  NAT ( »T  SEC*  ) 

CALL  TEXTGt I  DEV, I  TEXT, 24, 22, 75, 1,2,  IER) 

I F ( IER.NE.O) OUTPUT ( 10 1 ) IER, 'PTI ME  ' 

ENCODE (  96 ,203, ITEXT  )N ELK 
FCFNAT( 'STATISTICS  ON  BLOCK', 13) 

CALL  TEXTQ( IDEV, ITEXT ,24,  24 , 3 2 ,  1 , 2 , 1 ER ) 

IF (IER.NE.O  OUTPLT (101) I 5R,*TE1* 

ENCODE ( 96 ,204, ITEXT ) MEAN, SDEV 
FCF  NAT ( 'MEAN  VALUE= • , F8.3 , 4X  ,  ' S C EV=  •  , F8 
CALL  TEXTC( I  DEV, I  TEXT  ,24,2  5,24, 1 ,2  ,  IER) 

IF  (IER.NE.O  OUTPLT ( 101  )  IER,  'TE2' 


2  CO 


30 


3C4 


F8.3  ) 


lF(lEP.NE.G)  OUTPUT (101 _ 

ENCCCE(9c, 206, ITEXT)  MI N »  M  A  X 
206  FCFNAT('MIN  PERIOC  IN  NICSCSEC 
14 X  ,  'MAX  PEP IOD= ' , F8 .3  ) 


-  i 


F8.3 
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,  IVAR (6  I  , 


non  nno  or>o  oor> 


20  £ 

2  1C 

209 

42 

40 

41 
1 2  C 

215 

2  1  £ 

99E 


cm  TEXTC(IDEV»ITEX7 ,24,26,24, 1,2, IER) 

IF  < IER.NE.O)  OUTPUT  (101  HER,  'T52  • 

ENCCCS(  96, 308, ITEXT) 

FC  FMAT (  '  SET  MIN  £  MAX  VALUE  CF  Y  AXIS  IN  MICRSEC') 
CALL  7 EXTG ( I  DEV, IT  EXT, 24, 28, 24,  1,2,1 ER) 

IF(IER.NE.O)  OUTPUT  ( 101  )I£R,,TE4' 

ENCCDE(96, 210, ITEXT) 

FCFMA7(  ’TYPE  A  MIN, A  COMMA,  ANC  MAX  IN  AN  F  FORMAT  *  ) 
CALL  T£XTG( IDEV, ITEXT, 24, 29 .24, 1,2,  IER) 

I F ( IER.NE.O) OUTPUT (101  HER, JPTYPE' 

ENCCCE(96, 209, ITEXT) 

FC  F  MAT (  ' Tu  USE  STATISTICAL  MIN  ANC  MAX  AS  Y  AXIS 
IN  IN  ANC  MAX  TYPE  JUST  A  C/R') 

CALL  T5>T0(IDEV, ITEXT ,24, 22 ,24 , 1  ,2,IER) 


IFUeR.NE.C)  OUTPUT  (101)'  I  £  R ,  '  T  I 


24, 

i 


CALL  TEXTR( I  DEV,  77  777  7  776,1 ,20, 2 
I  F  (  I  ER • N  E • 0 ) OUTPUT ( 10 1 ) I ER , ' TXTR 
I  F  (MGC( ITDIR(9) ,8).EQ.O)  GO  TO  40 
CALL  TEXTI(IDEV,RESP,4,0, 9, IER ) 

1 F (PESP (1) .N6.0H  )  GC  TO  41 

CMIMMIN 

CM A>=MA> 

CCNTINLE 

CECCCEt 16, 120, RES P  )CMIN,DMAX 
FCFMAT(2F8.0) 

I F ( CM  IN .CE.DMAX)  GC  TO  42 
CM  AX2  =  C M  AX-CNIN 
ENCCCE( 96 ,215, ITEXT)  DM  IN 
FCPNAT<  F8.2> 

CALL  7EXT0( I  DEV, IT  EXT, 24, 20, 20, 1,2, 
IF( IER. NS.O)OUTPUT( 101) IER, • LAST  • 
ENCCCE ( 96,216, 1  TEXT)  DMA  X 
FCFMAT(F8.2) 

CALL  7E  >70 ( IDEV , ITEXT ,24, 5, 

IF (lES.NE.C)GUTPUT(lOl)  IER,  • 

K  K  K  =  1 
CCNTINLE 


1,2, IER ) 


2C, 1,2, IER) 
LAST  +  1 » 


ELFFERIN  4096  POINTS  CF  ARRAY  CCUNT  IMG  INCATA 
I  C  E  V=  1 

CALL  eUFFER IN (1,1, INDATA, NS  AMP,  INC) 

4  I  F  (  INC.  EC  •  1 )  GO  TG  4 

I  F  (  (N8LK.LE.il ) .CR. (NELK.GE. 14)  )  GC  TO  301 
CC  27  L 1  =  1 ,4081 , 16 

WFI7E(6,106) ( INDATA ( 12 ),L2»L1,L 1+15) 

1C6  FC F MAT ( IX, 16(16, 2X  )  ) 

27  CCNTINUE 
SCI  CCNTINUE 

MAKE  INCATA  A  REAL  ARRAY 

CC  10  1=1, N SAMP 

CATA( I ) =FLCAT (INOATA(I  )  )/4. 

1C  CCNTINUE 

IF ( (KKK.EC.3 ) .AND. (NEWI.Nc.L6LK ) )LREC=N£WI-2*N$AMF 
NCFMALIZE  (0  TO  1)  TC  CISPLAY 
CC  12  I =1 , LREC 

C  AT  A (I )  =  (CATA( I )-CMIN ) /CMAX2 
12  CCNTINUE 

FCFM  7NE  ARRAY  INDATA  WITH  GRAPHICS  WCRCS 


L  =  1 

996  CCNTINUE 

CALL  CC-IN  IT  ( IDEV,  IGCIR  ,5,  IER) 

IF(IER.ME.O)  OUTPUT (101)  IER,'D1' 

T  (1  )=TIME 
CC  19  I  *  1 , N D I S 

TIME  =  TIME+(( (CATA(L)*GMAX2) +DM IN  )  / 1000000.  ) 


88 


non 


I F ( I. EQ • ( NC I  S/2 )  ) T ( 2 ) =T IME 
IF(  I.EQ.NDIS)  TO  1  =T IMS 
IS  CCMINLE 

ENCODE ( S6 ,320, 1  TEXT )T(1),T(2),T(3) 

320  FCFf'AT<2  8X,F6.3,12X,F6.3,12X,F6.2) 

CALL  TEXT0<IDEV,ITEXT»24,22»l,l,i*IER) 

I  F (  IER. NE.O) OUTPUT (101) IER, 'TTIVE* 

INCATA(  1)=IFEAD(0, 10) 
1NCATA(2)=IPACK(-.5,.5,Q) 

INCATA(  3  )  =  IPACK(  (  1 .  /D  I  V  )-.  5  »CAT  A  (0,0) 

CC  6  I  =  2  »  N  C I S 

INCATA(  I+2)  =  IPACK( ( I /O I  V )-. 5 , DAT A ( L  +  I-l  )  ,1  ) 
6  CONTINUE 
NlfC  FD  =  N  C  IS  +  2 

CPAW  CN  SCREEN 

CALL  GRAPHG( IDEV, IAXIS ,9,1, IER ) 

I  F  (  1ER.N£.0)GUTPUT <101)IER, *G1» 

CALL  GRAFHC ( ICEV,  IN  CAT A .N WORC , 2, 1 EF ) 

IF  ( lER.NE.O)QUTPUT(lOl) IER, 'G2' 

N  =  N  4NC I S 

EC  CALL  FNS  (  ICEV  ,  I  Sin  ,  I  ER  ) 

IF(IER.NE.Q)  OUTPUT ( 101 )  IER, 'FI' 

IF  (LLS( ISU,15).LT.C )  GC  TO  18 
IF<LLS(ISW,16).LT*0»  GO  TC  995 
GC  TO  80 

is  l  =  l+nl;is 

I  F  {  L.  I E  •  ( N SAMP-NO I S+ 1 ) )  GO  TC  996 

KI<K  =  KKK-4l 

IF<KKK.GT.2)N8LK=NELK+1;GG  Tu  999 
GC  TO  SSS 
9 S  5  S TC  F 

Er*c 

-LOAD  XF,PAP 
>C  M  A 
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